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ABSTRACT 

We use data on the high-redshift evolution of the size distribution and luminosity function of galaxies 
to constrain the relationship between their star formation efficiency and starburst lifetime. Based on 
the derived scaling relations, we predict the angular sizes and average surface brightnesses of faint 
galaxies that will be discovered with JWST. We find that JWST will be able to resolve galaxies at 
the magnitude limit toab < 31 out to a redshift of z ~ 14. The next generation of large ground-based 
telescopes will resolve all galaxies discovered with JWST, provided they are sufficiently clumpy to 
enable detection above the bright thermal sky. We combine our constraints with simple models for 
self regulation of star formation, and show that feedback from supernovae at redshifts z > 3 is likely 
mediated through momentum transfer, with the starburst timescale set by the lifetime of the massive 
stars rather than the dynamical time in the host galactic disk. 

Subject headings: galaxies: evolution, formation, high-redshift — cosmology: theory 



1. INTRODUCTION 

The luminosity function of Lyman-break galaxy candi- 
dates discovered at z > 6 in the Hubble Ultra-Deep Field 
is described by a Schechter function with a characteris- 
tic luminosity that decreases towards higher redshift as 
expected from the dark matter halo mass-function (e.g. 
Munoz & Loeb 2010). At low luminosities, the lumi- 
nosity function is fit as a power-law with a faint end 
slope of a ~ —1.8 that is roughly independent of redshift 
(McLure et al. 2009; Bouwens et al. 2010; Yan et al. 
2010). The shape and overall density of the luminosity 
function is the primary observable that must be repro- 
duced by any successful model of galaxy formation (e.g. 
Raicevic, Theuns & Lacey 2010, Salvaterra, Ferrara k, 
Dayal 2010; Trenti et al. 2010). 

In a recent study, Oesch et al. (2009) have measured 
the redshift evolution of the scale-length of galactic disks 
(see also Ferguson et al. 2004; Bouwens et al. 2004). 
Among a sample of galaxies with constant luminosity, 
the half-light galaxy radius was found to scale as, 



-Rgal OC (1 + z)~ 



(1) 



with m = 1.12 ± 0.17. This evolution is consistent with 
the inverse relation between virial radius and redshift 
[i? ga i oc (1 + z)~ ], that is expected for dark-matter ha- 
los assuming a constant halo mass to luminosity ratio 
(Ferguson et al. 2004; Oesch et al. 2009). 

In this Letter we model the evolution of galactic disk 
radii for different parameterised models of self-regulated 
star formation. We parameterise our models in such a 
way that the constraints based on the observed size of 
galaxies are orthogonal to those derived from the slope 
of the galaxy luminosity function. Based on an empiri- 
cal determination of parameters of this model, we predict 
the expected size distribution of galaxies in future sur- 
veys with JWST, and show that existing observations al- 
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ready limit the physical processes that govern the global 
properties of star formation at high redshifts. Where 
required, we adopt the standard set of cosmological pa- 
rameters (Komatsu et al. 2009), with density parameters 
n h = 0.044, n m = 0.24 and fi A = 0.76 for the matter, 
baryon, and dark energy fractional density, respectively, 
and h — 0.72, for the dimensionless Hubble constant. 

2. MODEL 

To model high redshift galaxy properties we begin with 
an expression for the UV luminosity of a galaxy within 
a dark matter halo of mass Mh a i , 



L oc /*M ha i 



mm(t SB ,0 
tsB t± 



(2) 



where f+ is the fraction of baryons turned into stars, 
tsB is the lifetime of the starburst, and t+ is the aver- 
age lifetime of the massive stars (> 8M ) that affect the 
feedback through supernova (SN) explosions and domi- 
nate the UV luminosity output L. The lifetime of main- 
sequence stars, t ms = lO lo yr(M star /M ) -2 - 5 (Hansen, 
Kawaler & Steven 1994), implies an average delay time 
for supernova feedback of t* ~ 7 x 10 6 yr for an ini- 
tial mass function of massive stars dN/ dM stax oc M s "jT ar 35 
(Scalo 1986). The value of t$B is thought to be related 
to the dynamical time (idyn) of the galactic disk (Ken- 
nicutt 1998), which scales as the age of the Universe tn, 
namely t dyn - 3 x lO" 3 ^ « 3 x 10 6 yr[(l + z)/7]~ 3 / 2 . 
Equation (2) implies that the luminosity could become 
independent of the starburst lifetime at sufficiently high 
redshifts z > 3 for which tdyn < t*- The mass-to- light 
ratio is governed by /* and £sb, each of which can de- 
pend on both Mh a i and z. We therefore parameterise a 
general class of models for high redshift galaxy formation 
using the ratio, 



(3) 



Physically, this quantity is proportional to the inverse 
of the mass-to-light ratio. This parameterisation can be 
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Fig. 1. — Left Panel: Contours of parameter combinations (a, 6) that give the best fit observed value m = 1.12 (central solid line), as 
well as the ±l-sigma (outer solid lines) and ±2-sigma combinations (dashed lines). Central Panel: Contours of parameter combinations 
(a, 6) that give the best fit observed value a' = 0.8 (central solid line), as well as the ±l-sigma (outer solid lines) and ±2-sigma (dashed 
lines) combinations. Right Panel: The combined constraint from the observables m and a, derived from the product of likelihoods 
C m = exp[-(m- 1.12) 2 /2(0.17) 2 ] and C a = exp[-(a' - 0.8) 2 /2(0.08) 2 ]. The contours represent the 32%, 11% and 4.5% likelihoods, which 
for a Gaussian distribution contain 68%, 90% and 95.4% of the probability, corresponding to the 1-sigma, 1.5-sigma and 2-sigma ranges 
for parameter pairs. The model assumes g = 0. For comparison, the solid and open circles show the example where there is no feedback, 
but the starburst timescale is proportional to the dynamical time of the galaxy (a = 0, 6 = 3/2) or the average lifetime of massive stars 
(a = 0, b = 0), respectively. Similarly, the solid and open diamond points show the example of the model in Wyithe & Loeb (2003) that 
proposes star-formation is limited by the production of the binding energy of the galactic gas in the form of SN-driven winds with tgB > t* 
(a = 2/3, b = 5/2) and tgg < t* (i = 2/3, b = 1). Finally, the solid and open triangle points show the example of a model where star 
formation is limited by SN-driven winds which deposit momentum rather than energy into the galactic gas, with tgB > t* (a = 1/3, 6 = 2) 
and t SB < U (a = 1/3, 6 = 1/2). 



used to describe a range of physical models that predict 
the unknown variation of the star-formation efficiency 
and starburst lifetime with redshift and halo mass. Each 
particular model of the self-regulation of high redshift 
star formation will yield distinct values for the parame- 
ters a and b. 

In order to compare the resulting grid of models with 
the observed evolution in galaxy size we require two fur- 
ther ingredients, namely the virial relation (simplified 
form valid at high redshift), 



halo 



cx i? 3 ir (l + zf 



(4) 



and the parameterisation, 



-Rgal OC i?vir(l 



to describe the relation between virial radius i? v „- and the 
galactic disk scale radius i? ga i- In the latter expression 
we expect a parameterisation with g = if the disk size 3 
is set by the adiabatic model of Mo, Mao & White (1998). 
The above set of relations can be used to find the pre- 
dicted evolution of the galaxy radius with redshift at a 
fixed luminosity, yielding 



m = 1 — g 



1 b 



3a+ 1 



(6) 



which can be compared with the observed value of m = 
1.12 ± 0.17 (Oesch et al. 2009). 

It is also possible to use the observed slope of the lu- 
minosity function to constrain the model, thus break- 

3 If the gas disk becomes stable to fragmentation at a radius 
beyond which there is a significant fraction of gas by mass, then the 
half mass radius of the stellar disk may not equal the scale radius of 
the gas-disk. However, we find that the disk becomes stable (based 
on Toomre Q criterion) only at 3-4 scale radii, and therefore that 
there is a very small fraction of mass at these large radii. Thus, we 
adopt the assumption g = throughout our analysis. 



ing the degeneracy between the parameters a and b that 
arises from application of equation (6). Modelling the 
galaxy luminosity function using the halo mass func- 
tion (dn/tiMhaio) and the simple parameterised model 
described in equations (2-3) above we find 



a ■ 



d log n d log n 



cHogL dlogMhaio 

1 d log n 
a + 1 d log M ha i ' 



<f l0gM ha io 



dlogi 



(7) 



In the mass range 10 8 M Q < M halo < 10 10 5 M Q (near 
or below the non- linear mass scale at z ~ 6), the 
mass function has the local power-law slope —2.05 > 
W d log n/d log Mhaio ^ —2.45. At low luminosities the 



Ihal 

high redshift galaxy luminosity function is fit as a power- 
law with faint end slope a — —1.8 (McLure et al. 2009; 
Bouwens et al. 2010; Yan et al. 2010). While the un- 
certainty in a is large at z ^ 7, it is well constrained in 
the redshift range 4 < z < 6. We estimate a 10% uncer- 
tainty in the value of (d log n/d log Mhaio), dominated by 
the uncertain mass of the host halos. We therefore define 
a' = a x (d log n/d log Mhaio ) ~ 1 , an d adopt the constraint 
a' = l/(a + 1) = 0.8 ± 0.1 based on equation (7). 

2.1. Parameter constraints 

The left hand panel of Figure 1 shows contours of pa- 
rameter combinations (a, b) that give the best fit ob- 
served value to = 1.12, as well as the ±l-sigma and 
±2-sigma combinations (equation 6 assuming g = 0). 
The central panel of Figure 1 shows contours of param- 
eter combinations (a, b) that give the best fit observed 
value a' = 0.8, as well as the ±l-sigma and ±2-sigma 
combinations (again assuming g — 0). The right hand 
panel shows the combined constraint from both observ- 
ables, with the contours representing the 32%, 11% and 
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Fig. 2. — The projected relation between galaxy size and redshift, plotted for four values of apparent magnitude. The Left and Central 
panels show the physical (H ga i) and apparent angular sizes (6 ga i) respectively. For comparison (thick grey lines), we also plot the indicative 
angular resolution A6* of telescopes with diameters corresponding to HST (D tc \ = 2.5 m), JWST (D tc \ = 6.5m) and ELT (D tc i = 30 m). 
The Right panel shows the relation between surface brightness (averaged within the scale radius) as a function of redshift. In each panel, 
the grey band around the case of ttiab = 29 mag shows the 68% range of uncertainty on the mean. 



4.5% likelihoods, which for a Gaussian distribution cor- 
respond the 1-sigma, 1.5-sigma and 2-sigma regions for 
parameters combinations. We find a — 0.3 ±0.15 and 
b = 0.5 ± 0.6 (68% range on individual parameters). 

3. THE SIZE OF HIGH REDSHIFT GALAXIES 

Our empirical constraints can be used to extrapolate 
galaxy size out to the higher redshifts and lower lumi- 
nosities that will be observed by the next generation of 
telescopes. Beginning with equations (2-3) we get, 

LocM^l + z) 6 , (8) 
which when combined with equations (4-5) implies, 
R ga i(xL^+^ (1 + z )-( 1 -9+3(tW) a L3(iW(l + z )- m 



and an extremely large telescope (ELT; D tc \ — 30 m) 

-l 



/ D ( z \ \ m AB 1 Q~ m AB f 1 + Z\ ' 

-*(b;Kj) 10 (its) < 9 > 

where Z?l is the luminosity distance, toab is the appar- 
ent galaxy magnitude, and mAB.o is the apparent mag- 
nitude of galaxies from which the normalization of the 
relation (Rq) is measured at redshift z a . The normalisa- 
tion of this relation is calibrated using observed galaxy 
sizes at z ~ 7—8 (Oesch et al. 2009). We average 
the prediction over the three independent points 4 having 
(z ,i?o,m AB ,o) = (7, 0.5 ±0.1, 27.8), (7, 0.75 ± 0.1, 26.6), 
and (8,0.4±0.1, 28.1), and over distributions to = 1.12± 
0.17 and a = 0.3 ± 0.15. The resulting relation is shown 
in Figure 2 as a function of redshift for four values of ap- 
parent magnitude. In the case of toab = 29 mag, the grey 
band shows the 68% range of uncertainty on the mean 
radius. The left and central panels show the physical 
and apparent angular sizes (# ga i), respectively. Galaxies 
of fixed apparent magnitude have smaller physical sizes 
(but larger angular sizes) at higher redshift. 

For comparison, we also plot in Figure 2, the indica- 
tive angular resolutions of telescopes with diameters cor- 
responding to HST (Ad = 2.5m), JWST (D tcl = 6.5m) 

4 Uncertainty in Rq is the standard error on the mean. 



(10) 



Here we have calculated the resolution at the observed 
wavelength of the rcdshiftcd Lya line A = 1216(1 + z)A. 
This simple comparison suggests that galaxies with an 
apparent magnitude of itiab = 28 mag (1 magnitude 
fainter than the HST WFC3/IR limit) have an angular 
size of 9 ga \ ~ 0.1", which is at the limit of HST resolu- 
tion (Oesch et al. 2009). Fainter galaxies at higher red- 
shifts cannot be resolved by HST. The larger aperture of 
JWST will allow the study of galaxy structure at higher 
redshift and fainter fluxes (Windhorst et al. 2008). For 
example, galaxies with toab = 29 mag will be resolved 
out to at least z > 15. Moreover, JWST will be able 
to resolve galaxies at the magnitude limits of toab < 30 
and toab < 31 out to z ~ 10 and z ~ 14, respectively. 
Galaxies at yet higher redshifts or fainter fluxes would 
appear point-like. Bright galaxies are known to become 
rarer toward high redshift (e.g. Bouwens et al. 2010; Yan 
et al. 2010), and the discovery of such very high redshift 
galaxies by JWST could be limited by its relatively small 
field of view (e.g. Munoz & Loeb 2010). However, our 
extrapolation suggests that an ELT would be able to re- 
solve all galaxies detectable with JWST. 

The above discussion refers only to angular resolution 
and neglects surface brightness sensitivity. Calculation 
of the details of surface brightness sensitivity (see, e.g. 
Windhorst et al. 2008) are beyond the scope of this Let- 
ter. However, we note that faint, high-redshift galax- 
ies are always below the surface brightness of the sky. 
Therefore, in order to measure the details of a high red- 
shift galaxy profile, an HST observation must achieve a 
very high signal-to-noise ratio on the zodiacal sky (of 
order 10 3 per pixel) so that the background would be re- 
liably subtracted (Hathi et al. 2008). The right panel of 
Figure 2 shows the average surface-brightness within the 
scale-radius of galaxies of various magnitudes as a func- 
tion of redshift. Galaxies with toab = 29 mag show ~ 
24.5 magnitudes per square arcsecond at z ~ 6, becoming 
fainter by only a modest amount towards z ~ 20. The 
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plotted curves should be compared to the space based 
value for the zodiacal sky in the K-band of ~ 22.5 mag- 
nitudes per square arcsecond (Windhorst et al. 2010), or 
the ground based K-band value of ~ 14 magnitudes per 
square arcsecond for the thermal sky. The difference in 
sky brightness between the ground and space implies that 
HST is equivalent in sensitivity (but not in resolution) to 
a 150 m ground-based telescope for the purpose of imag- 
ing resolved high redshift galaxies. Thus, a ground-based 
ELT will not be competitive with HST or JWST, since it 
must overcome the much brighter sky. As a result, even 
though all high redshift galaxies discovered with JWST 
could be resolved by an ELT, their surface brightness will 
render their extended structure undetectable. However, 
high redshift galaxies are observed to be very clumpy 
owing to the presence of star forming regions (Hathi et 
al. 2008), and an ELT wi\\ be more sensitive to these un- 
resolved clumps than to diffuse structure, owing to the 
higher resolution in addition to larger collecting area. In- 
deed, HST is comparable to only a 20 m ground based 
telescope with respect to point source sensitivity (but 
without the resolution of a 20 m aperture). Thus, if high 
redshift galaxies are sufficiently clumpy, then an ELT 
may be able to obtain close to the theoretical resolution 
shown in the central panel of Figure 2. 

4. STAR FORMATION MODEL CONSTRAINTS 

A range of simple models for the process governing star 
formation can be compared to the constraints on param- 
eters a and b. We consider three simple models for the 
possible feedback. The first is a model where there is 
no self regulation of star formation. The second and 
third models describe the evolution in cases where the 
star formation is limited by SN feedback, and the inter- 
action between the SN driven wind and the galactic gas 
conserves energy and momentum, respectively. In each of 
these three cases we consider scenarios where the lifetime 
of the starburst is proportional to the galaxy dynamical 
time and to the lifetime of massive stars, respectively. 
Altogether, we have 6 model predictions for m and a' 
with which to compare our constraints. 

4.1. Models without feedback 

We begin with feedback-free models in which the star 
formation efficiency /* cx M*/Mh a io = const. We first 
consider the case where the natural timescale for the 
starburst is proportional to the galaxy dynamical time, 
tsB oc ^dyn oc (1 + z)~ 3 / 2 , which implies 

min(t SB ,M A sc M h°aio( 1 + ^)° M o (1+z) 3/2 

yielding a — and b = 3/2. This parameter combination 
is shown by the solid circles in Figure 1 for comparison 
with the present constraints. We find that this simple 
model is inconsistent with both the observed evolution 
in i? ga i and the observed luminosity function slope at the 
2-sigma level. The combined constraint rules out this 
model at high confidence. 

However, at the high redshifts when the dynamical 
time is shorter than the average lifetime of massive stars, 
the same luminosity may be achieved with a lower star 



formation rate than in the model described above, 

/ ^f^ ocf ocM° alo (l + ,)°, (12) 

tSB t* t* 

yielding a = b = for galaxies in this case. Such a 
model (shown by the open circular symbols) is ruled out 
by the slope of the luminosity function at the 2-sigma 
level, but is consistent with the observed evolution in 
galaxy radius. The combined constraint rules out this 
model at high confidence. 

4.2. Models with feedback through energy conservation 

We next consider models including self-regulation of 
star formation by supernovae feedback. The model of 
Wyithe & Loeb (2003) assumes that star formation is 
limited by the transfer of energy from SN-driven winds 
when it is equal to the binding energy to the galactic gas 
(Dekel & Woo 2002). The model asserts that stars form 
with an efficiency /* out of the gas that collapses and 
cools within a dark matter halo and that a fraction Fsn 
of each supernova energy output, Esn, heats the galac- 
tic gas mechanically. The mechanical feedback will halt 
the star formation once the cumulative energy returned 
to the gas by supernovae equals its binding energy (as- 
suming negligible radiative losses for a sufficiently rapid 
starburst). Hence, in this model the limiting stellar mass 
is set by the condition 



EsnFsn 



E h 



1 fit 

2 a, 



-M ha i Q v c 



(13) 



where wsn is the mass in stars (in M©) per supernova 
explosion. Equations (4) and (13) imply that the total 
mass in stars, = (/*O b /f2 m )-Mh alo , scales as 

M*«A^(1 + *), (14) 

and hence the star formation efficiency scales as /* oc 

My'(l + z). Thus, smaller galaxies are less efficient at 
forming stars, but a galaxy of fixed mass is more efficient 
at forming stars at higher redshift. 

In a model where the starburst lifetime is proportional 
to the galaxy dynamical time, we find 



halov ~ ' oc M, 



dyn 



1 



-3/2 



2/3 (1 
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N5/2 



(15) 



This model can therefore be parameterised in terms of 
the combination a = 2/3, and b = 5/2 for galaxies, which 
is shown by the solid diamonds in Figure 1. This simple 
supernovae driven feedback model is inconsistent with 
the observed evolution in i? ga i and the luminosity func- 
tion slope at the 2-sigma level. The combined constraint 
rules out this model at high confidence. 

However, at high redshifts z > 3 when the disk dy- 
namical time is shorter than the lifetime of a massive 
star t+ ~ 10 7 yr, we note that SN feedback will be less ef- 
ficient, and the star formation efficiency could exceed the 
self regulated value of /* described above. In particular, 



«M*£(1 + *), 



(16) 



yielding a — 2/3, and b = 1 for galaxies in this case (open 
diamond symbols). Such a model is consistent with the 



■5 



observed evolution of galaxy size, but inconsistent with 
the slope of the luminosity function at the 2-sigma level. 

4.3. Models with feedback through momentum 
conservation 

Finally, we consider a model where the SN-driven 
winds conserve momentum in their interaction with the 
galactic gas rather than energy. In this case the limiting 
stellar mass is set by the condition 



w SN c 



SN 



1 fit 



-M ha \ Q v c . 



(17) 



In a model where the starburst timescale is proportional 
to the galaxy dynamical time, we find 



kdyn 



(18) 



This model (solid triangles), represented by a = 1/3 and 
b = 2, is consistent with the constraints from a, but 
ruled out at the 2-sigma level by the constraints from 
the evolution of galaxy radius. If instead the lifetime of 
massive stars sets the starburst timescale, we find 



/* 



a/2 



(19) 



which is represented by a = 1/3 and b = 1/2, and is 
plotted as the open triangles in Figure 1. This model is 
consistent with both the constraints from m and a' . 

Our results imply that SN feedback in high redshift 
galaxies occurs through the transfer of momentum be- 
tween the SN-driven wind and the galactic gas, and that 



the starburst timescale is set by the average lifetime of 
massive stars rather than the dynamical time of the host 
galactic disk. 

5. DISCUSSION 

We have used the observed redshift evolution of disk 
sizes, and luminosity function of galaxies to constrain the 
relationship between star formation efficiency and star- 
burst lifetime. We find that supernova feedback in high 
redshift galaxies is likely mediated through momentum 
transfer with the starburst timescale dictated by the av- 
erage lifetime of the massive stars, t*. The latter result 
follows naturally from the fact that the dynamical time 
of galactic disks becomes shorter than t+ ~ 10 7 yr at 
redshifts z > 6. 

We extrapolated the derived scaling relations to pre- 
dict the angular sizes of galaxies at higher redshifts 
and fainter fluxes than currently observed. We have 
found that JWST will be able to resolve galaxies with 
TTiAB < 30 and ttiab < 31 only out to redshifts of z ~ 10 
and z ~ 14, respectively. However, if high redshift galax- 
ies are sufficiently clumpy, so that unresolved star form- 
ing regions can be detected above the bright thermal sky, 
then the next generation of ground-based extremely large 
telescopes will be able to resolve structure in all galaxies 
at all redshifts detectable by JWST. 
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